Revised Phase Diagram of the Gross-Neveu Model 
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We confirm earlier hints that the conventional phase diagram of the discrete chiral Gross-Neveu 
model in the large N limit is deficient at non-zero chemical potential. We present the corrected phase 
diagram constructed in mean field theory. It has three different phases, including a kink-antikink 
crystal phase. All transitions are second order. The driving mechanism for the new structure of 
baryonic matter in the Gross-Neveu model is an Overhauser type instability with gap formation at 
the Fermi surface. 



I. INTRODUCTION 

The simplest variant of the Gross-Neveu (GN) model is 
a 1+1 dimensional relativistic field theory with N species 
of fermions interacting via a quartic self-interaction [1] , 



C = ^d^ + \g 2 {W) 



(1.1) 



where we suppress the flavour indices (ipip = J2i i^i^i 
etc.). This model possesses a discrete chiral symmetry 
(ip — > 7 5- 0) which forbids a bare mass term. In the 
large N limit, it can be solved exactly. It exhibits a 
surprisingly large number of phenomena of interest to 
nuclear and particle physics (for a recent review, see Ref. 
[2]). Asymptotic freedom, dimensional transmutation, 
spontaneous breakdown of chiral symmetry, dynamical 
mass generation, quark-antiquark and multiquark bound 
states are its most prominent features. Moreover, the 
GN model is believed to have a non-trivial phase di- 
agram. The discrete chiral symmetry gets restored at 
high temperatures and/or densities. According to com- 
mon lore, a critical line which separates the chirally sym- 
metric phase (massless fermions) from the broken phase 
(massive fermions) comprises both first and second or- 
der lines separated by a tricritical point [3]. Models for 
which one can construct the renormalized phase diagram 
analytically are extremely scarce and worth studying on 
theoretical grounds, even if they are far from being re- 
alistic. In the present case it is essential to define the 
model by letting N — ► oo before taking the thermody- 
namic limit. In this way one can avoid conflict with no- 
go theorems characteristic for low dimensional systems 
[4, 5]. In particular, spontaneous breaking of continu- 
ous symmetries becomes viable in 1+1 dimensions [6]. 
In spite of all the obvious reservations we feel that one 
can learn more about relativistic matter from Eq. (1.1) 
than one had any right to hope on the basis of such an 
incredibly simple model. 

Given the large number of papers dealing with the GN 
model, one may wonder whether there are any open ques- 
tions left. This is actually the case. In the present work, 
we shall reconsider the model at finite temperature and 



chemical potential. As we have pointed out recently, the 
alleged phase diagram of the GN model is likely to be 
flawed [2]. It does not do justice to the role played by 
baryons in the structure of baryonic matter through ef- 
fects which are not 1/N suppressed. Formally, the prob- 
lem with previous investigations lies in the tacit assump- 
tion of unbroken translational invariance. The conden- 
sate (i/jtp) then gives rise merely to an x-independent 
mass term. Following similar studies of the GN model 
with continuous chiral symmetry [7, 8], we have found 
evidence for a crystal phase with lower energy than the 
standard massive Fermi gas at zero temperature — a 
kind of kink-antikink crystal [9]. This result was based 
on a variational calculation using a specific ansatz. In 
the present work, we report on a complete Hartree-Fock 
calculation at finite T and /i. This enables us to round 
off our series of investigations by presenting the revised 
phase diagram of the discrete chiral GN model in "full 
glory" . The basic theoretical tool we have employed here 
is the relativistic Hartree-Fock method, appropriate for 
the leading order term of the large N expansion. The 
same method has been used in Refs. [10], [7] and [9] to 
which we refer the reader for further background mate- 
rial. 

This paper is organized as follows: In Sect. II, we ex- 
plain how we have solved the Dirac-Hartree-Fock equa- 
tion using a combination of numerical and analytical 
methods. Sect. Ill is dedicated to cold baryonic matter 
(T = 0). In Sect. IV, we collect the necessary formu- 
lae for finite temperature and chemical potential calcu- 
lations. Sect. V contains the main result of this work, 
namely the revised phase diagram of the GN model. In 
Sect. VI, we summarize our findings, compare them with 
the continuous chiral GN model and put them into per- 
spective. 



II. SPECTRUM OF THE SINGLE PARTICLE 
DIRAC HAMILTONIAN 



We first discuss how to solve the Dirac-Hartree-Fock 
equation. In the GN model with (-0?/;) 2 -interaction, a 



2 



mean field approach can only generate a (local) Lorentz 
scalar potential. Therefore we have to diagonalize a sin- 
gle particle Dirac Hamiltonian of the generic form 

H = H + V = 1 5 ^+ 1 °S(x) . (2.1) 
1 ox 

We choose the following representation of the 7-matrices, 



7 = -ci 



7 1 = i°3 , 7 5 = -02 



(2.2) 



For spatially constant S(x), H reduces to a free, massive 
Dirac Hamiltonian and our task would be rather trivial. 
This would evidently lead back to the standard results 
for the GN model. Instead, we admit the possibility that 
translational invariance breaks down and allow for spa- 
tially periodic S(x) with (yet to be determined) period 
a, 



S{x + a) = S(x) . 



(2.3) 



We first diagonalize the free Hamiltonian Ho. This step 
serves both to generate a basis for the full, numerical 
diagonalization of H and as a starting point for pertur- 
bation theory. Enclosing the system in a box of length 
L = Na and imposing antiperiodic boundary conditions, 
we find (77 denotes the sign of the energy) 

H \ri,n)=ri\k n \\ri,n) (r? = ±1) , (2.4) 

with the normalized spinor wave functions 

(x\n,n) = -^=( . 1 ^ e ik " x , k n = HE( n +l/2). 

(2.5) 

Matrix elements of the potential V = ■y°S(x) in the un- 
perturbed basis can be expressed via the Fourier compo- 
nents of S(x) 



S(x) 



'See 



i2ir£x/a 



(S- t = s* e ) (2.6) 



as follows, 



(7]',n'\V\r],n) = - [risgn(fc n ) - r/sgn(fc„,)] S n 



(2.7) 



Eqs. (2.4) and (2.7) are all that is needed for the numer- 
ical diagonalization of H. Due to the spatial periodicity 
of S(x), the Hamiltonian matrix assumes a block diag- 
onal form (conservation of Bloch momentum). Further 
simplifications occur due to additional conservation laws 
for certain types of potential. Most of our calculations 
are based on the additional assumptions 



S(-x) - -S(x) 
S(x + a/2) = -S(x) 



(2.8) 



Such a scalar potential consists of a regular succession 
of positive and negative bumps of identical, symmet- 
ric shape. For this choice, only odd I appear, the Sg 



are purely imaginary and the Hamiltonian matrix is real 
symmetric, cf. Eq. (2.7). In some cases, we have em- 
ployed more general shapes of the potential which then 
lead to larger, complex matrices. It turns out though that 
the self-consistent Hartree-Fock potential belongs to the 
more restricted class of Eq. (2.8). 

A characteristic difficulty of the relativistic approach 
is the fact that one has to determine the Dirac sea self- 
consistently, not just the valence levels. In observables 
like the energy density or the baryon density infinities 
are then unavoidable. They require regularization and 
rcnormalization and hence a certain amount of analyt- 
ical work. Fortunately, due to asymptotic freedom of 
the GN model, fermion states with large momenta are 
only weakly affected by the mean field. Thus all UV- 
divergencies can be handled pcrturbatively. Perturbation 
theory has turned out to be very useful to the present 
study for other reasons as well. It provides more ana- 
lytical insight into various aspects of the phase diagram, 
in particular the vicinity of the chirally restored phase 
where S(x) is small. Therefore, we briefly outline how 
to apply perturbation theory to the Hamiltonian (2.1), 
restricting ourselves to potentials satisfying Eqs. (2.8). 

Naive application of 2nd order perturbation theory in 
V would yield the following energies, 



E v , n ~ r/sgn(fc„) [k n +^2 



2{k n - qg) 



with 



qg = — 



(2.9) 



(2.10) 



Eq. (2.9) displays simple poles which signal the break- 
down of naive perturbation theory in the vicinity of the 
gaps. As is well known from solid state physics, this 
problem can be cured by resorting to "almost degener- 
ate perturbation theory" (ADPT) in the region of the 
gaps [11]. More precisely, if Sg 7^ and we focus on the 
momentum region around qe, we have to diagonalize H 
exactly in two-dimensional subspaces of states whose mo- 
menta differ by a reciprocal lattice vector (2qg) but who 
have similar energy. For t — 0, these are states with op- 
posite 77 and the same (small) momentum k n ; for £ ^ 0, 
states with the same 77 but momenta k n ± qg, where k n 
is again small. The result can be stated concisely as fol- 
lows: Near qg, replace the term which blows up in Eq. 
(2.9) as follows, 

l^l 2 



— > sga(k n -q e )y/ (k n - qg) 2 + \Sg\ 2 -(k n -qg). 
2{k n - qg) 

(2.11) 

If one seeks a reasonable "global" approximation to the 
spectrum of H , one should apply Eq. (2.11) in the region 
between the midpoints of two gaps, i.e., 



k n e Ig 



qg + qg-i qg+i + qg 



(2.12) 



3 



Hence the perturbative eigenvalues are best denned 
piecewise as 

E $n ~ VSgn(k n ) (q e + sgn(/c„ - q t ) y/(k n - q e ) 2 + \S e \ 2 



\Si\ 



2(k n - qj) 



{k n e h) ■ 



(2.13) 



As is well known, perturbatively each Se ^ induces 
exactly one gap at k n w qt with width 2|S^|. 

Our last remark concerns the standard solution of the 
GN model which corresponds in our notation to So =/= 
only. In this particular case, ADPT yields the exact spec- 
trum, since the interaction only mixes the unperturbed 
states pairwise. 



III. BARYONIC MATTER AT ZERO 
TEMPERATURE 

We first present our results for the ground state of the 
GN model at finite baryon density (T = 0). We have 
to compute the Hartree-Fock energy density for a scalar 
potential S(x), 

For standard baryonic matter, the sum over occupied 
states includes all negative energy states (the Dirac sea) 
and as many positive energy states as required by the pre- 
scribed baryon density. Alternatively, one can consider 
"antimatter" where an appropriate number of negative 
energy states is unoccupied. We adopt here the latter, 
more economical way of evaluating the energy density. 
The 2nd term in Eq. (3.1) corrects for double counting 
of the interaction energy in the sum over single particle 
energies. Since the self-consistent potential must be of 
the form of a local, scalar potential in the GN model, it 
is sufficient to minimize the energy (3.1) with respect to 
S(x) in order to find the full Hartree-Fock solution. 

Eq. (3.1) is not yet in a form suitable for numerical 
calculations. The sum is badly UV divergent and the 
bare coupling constant appears in the 2nd term. Let us 
isolate the UV divergent terms in the sum as follows: The 
sum over negative energy states with short wavelengths 
is evaluated using perturbation theory. Provided we stay 
below the energy of the last gap, we can simply use Eq. 
(2.9), 



2 ^ 

L ^ E ~ hn 

n—n 



(3.2) 




2tt 



2tt 



A 



(k = k n + n/L). The quadratic divergence is irrelevant 
and can be eliminated by subtraction. The logarithmic 



divergence is cancelled by a corresponding divergence in 
the double-counting correction as can be shown with the 
help of the gap equation at zero density [10] (we choose 
units such that the vacuum fermion mass m is equal to 
1), 



Ng 2 



= 2 



dfc 



2tt TP 



= = - ln(2A) 

1 7T 



(3.3) 



The final expression for the regularized and renormalizcd 
ground state energy of a system with Fermi momentum 
labelled by n-p is then 



n=n F l>0 

(3.4) 

Before turning to the results, we briefly discuss two lim- 
iting cases where some analytical insight can be reached: 
The high density and the low density limit. 

In the high density limit, perturbation theory becomes 
applicable by virtue of asymptotic freedom. The only 
self-consistent perturbative solution we could find is one 
where S±i ^ only and all orbits below the first gap are 
filled (fcp = qi = n/a). Using Eq. (2.13) and the real 
variable 



Si = iSi , 

the self-consistent potential reads 

S(x) = 2Si sin(27rx/a) 



(3.5) 



(3.6) 



Performing the renormalization along the lines described 
above leaves us with 



£hf 



kp 

2-K ^ 47T 



31 



ln(l6k F Sl) - 1 



(3.7) 



where the first term is the result for a free Fermi gas of 
massless quarks. Minimization with respect to Si yields 
the equation 



Si ln(4fc/Si) =0 
with the two solutions 

1 



Si = 0, 



4fc F 



The non-trivial solution is lower in energy by 

1 



_ (£hf\ 



647Tfc2 



(3.8) 



(3.9) 



(3.10) 



Evidently at T = the discrete chiral symmetry does not 
get restored at any finite density. 

In the low density limit, S(x) is expected to be given 
by well separated, equidistant kinks and antikinks. We 
therefore approximate S(x) by tanhx for < x < A 
and 1 for A < x < a/4, then continue symmetrically. 
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Evaluating the Fourier coefficients Sg for such a potential, 
we obtain (in the limit A — > oo) 



Sg 



2k F 



sinh(7r^/cF) 



(3.11) 



The Fourier amplitudes satisfy some kind of scaling law, 
ISg being a "universal function" of Ikp. Roughly, Sg ~ 
l/£ and the convergence of the Fourier expansion of S(x) 
becomes quite slow. 

We now turn to the numerical results at T = 0. Wc 
have carried out the above sketched diagonalization pro- 
cedure. The renormalized energy density was then min- 
imized with respect to the period a and the Fourier co- 
efficients Si of the periodic scalar potential. We typi- 
cally work in a box of length L = 1000 (in units where 
the dynamical fermion mass in the vacuum is 1), tak- 
ing into account Fourier components in Eq. (2.6) up to 
I = 23 if necessary. As far as the period a is concerned, 
we confirmed the results of Ref. [9] : The system invari- 
ably adjusts its period such that the Fermi surface coin- 
cides with the bottom of the first gap. This is of course 
reminiscent of the Overhauser effect in condensed mat- 
ter physics [12] (for recent discussions in the context of 
QCD, see Refs. [13-16]). The spatially averaged baryon 
density per flavour is then given by 



P_ 

N 



k-p 

7T 



1 



(3.12) 



In Fig. (1) we plot the most important Fourier compo- 
nent (Si) as a function of k-p and compare it with the two 
analytical limits discussed above. The numerically found 
Fourier component interpolates nicely between the low 
and high density limits, a good test of the diagonalization 
and minimization procedure. Inclusion of higher Fourier 
components (I = 3,5, ...) is necessary at lower densities. 
After minimization of the energy the corresponding Sg 
are practically indistinguishable from the analytical re- 
sult, Eq. (3.11), so that we refrain from showing any plot. 
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FIG. 3: Hartree-Fock energy density compared with analyti- 
cal approximation obtained by asymptotic matching. 



A representative selection of the self-consistent Hartree- 
Fock potentials S(x) is shown in Fig. (2) for a range 
of densities. This figure illustrates how S(x) interpolates 
between well separated kinks and antikinks (<~ tanh x) at 
low density and the perturbative potential (<~ sin 2kpx) 
at high density. 

The energy density obtained from our numerical cal- 
culation is displayed in Fig. (3). The dashed curve is 
an analytical approximation obtained by matching the 
asymptotic behaviour at large and small kp as follows: 
For low density, we have widely separated baryons of 
mass 2N/n 7 therefore 



1 2k F 
N ~ 4tt tt 2 



{k F -» 0) . (3.13) 
At high density, perturbation theory predicts [cf. Eqs. 
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FIG. 4: Ground state energy difference between kink-antikink FIG. 5: Mean baryon density versus chemical potential at 



crystal and translationally invariant solution. 



T = 0, compared to the translationally invariant case. 



(3.7) and (3.10)] 
£hf 



fcp 



2tt 



647rfc|, 



(fcp — » 00 ) 



(3.14) 



The two expressions (3.13) and (3.14) have been matched 
at the point where they agree (fcp = 0.55537). This sim- 
ple estimate seems to catch the essence of a lengthy nu- 
merical calculation amazingly well. 

In Fig. (4), we show the energy difference between our 
solution and the standard solution of the GN model based 
on unbroken translational invariance (the dashed line is 
the perturbative prediction). The kink-antikink crystal 
is energetically favoured at all densities. Our results are 
of course also lower than the more restricted variational 
calculation of the kink-antikink crystal in Ref. [9], the 
improvement being most significant at high densities. Fi- 
nally, we can derive the baryon chemical potential from 
the slope of the curve in Fig. (3), 



d£ 



p 



HF 



dp 



(3.15) 



In Fig. (5) we have plotted p/N vs. p. As expected 
on the basis of the low density limit, the curve starts 
at Mb/N — 2/tt. It then rises steeply and approaches 
an asymptote. The solid curve is the result of the con- 
ventional calculation which involves a Maxwell-type con- 
struction. Here, the density jumps discontinuously at 
p = l/y/2 where a first order phase transition takes place. 
Our calculation is strongly in favour of a 2nd order phase 
transition at p = 2/n (or p = 0). These results at T = 
already indicate that the commonly accepted phase dia- 
gram of the GN model is not tenable. 



IV. GRAND CANONICAL POTENTIAL 

In this section, we collect the formulae needed for de- 
termining the phase diagram and computing thermody- 



namic observables in the grand canonical framework (see 
e.g. [17]). Starting point for the full, non-perturbative 
calculation is the formal expression for the grand canon- 
ical potential density <J> ((3 = 1/T), 



N 



A £ In [(l + e"^-")) (l + e^+^) 



(4.1) 



This quantity has to be minimized with respect to S(x). 
From now on, single particle energies will be defined as 
positive energies, i.e., E n = E + \ tTl . The sums over pos- 
itive and negative energy states correspond to the two 
factors under the logarithm. The first term in Eq. (4.1) 
is easily recognized as the "effective potential" expected 
for independent fcrmions. The 2nd term is identical to 
the double counting correction at zero temperature. Eq. 
(4.1) is still symbolic due to UV divergencies and re- 
quires some modifications. We proceed exactly like in 
the T — case, dividing up the spectrum into two parts: 
For n — 0...n, the energies E n are obtained by diagonal- 
ization of the Dirac-Hartree-Fock Hamiltonian in a large 
but finite box. For n > n we use the perturbative sin- 
gle particle energies and replace the sum by an integral. 
Using the notation 



e(k) 



E 



15/ 



1 



k - q e k + q e 



(4.2) 



for the perturbative energies and k as defined after Eq. 
(3.2), we find 



^ A £ In [(l + c-^-^) (l + e ^"+^) 

P n=0 

2 r°° dk 



N pL 

2 r 



In [(l + e-^*)-")) (l + e-^W+rt) 



6 



pk 

7T 



(4.3) 



t>o 



Note carefully the signs in the various exponents. They 
arise because one has to pull out the sum over single 
particle energies for k > k in order to perform the renor- 
malization, cf. Sect. III. As usual we have dropped the 
irrelevant divergent terms 



A2 
'2tt 



Ap 

TT 



(4.4) 



The quadratically divergent term is already familiar from 
T = 0, the linearly divergent term reflects the infinite 
baryon density of the Dirac sea. From Eq. (4.3), all other 
thermodynamic quantities (notably pressure P, baryon 
density p, entropy density s, and energy density u) follow 
in the standard way, 



2 — * 

1 80 



dp 

Ts-P + pp . (4.5) 



In the chirally restored phase (S(x) = 0), everything can 
be trivially evaluated in closed form, 



_ = It 2 + — - — — 

2tt 4tt ' N 
27r + 47r' AT 



A 
u 
A 



-T 2 + 



7T 



T . (4.6) 



V. REVISED PHASE DIAGRAM 

According to common lore, the GN model has two 
phases [3]: A chirally broken phase with a dynamical 
fermion mass m ^ at low (p,T), and a chirally sym- 
metric phase with massless fermions at high p and/or 
T (chiral symmetry referring to the discrete symmetry 
ip — > 7 5 i/>). These two phases are supposedly separated 
by a 2nd order line going from (p = 0,T = T c ) to 
(p = pt,T = T t ) and a 1st order line from the lat- 
ter point to (p = 1/V2, T = 0). The critical tem- 
perature at zero chemical potential has the value T c — 
c c /tt = 0.56693, whereas the tricritical point is located 
at p t = 0.60822, T t = 0.31833. A mixed phase also ap- 
pears which can be mapped out in other types of phase 
diagrams like the (u,p)- or (P, l/p)-plots, cf. Rcf. [18]. 

At variance with these results we have clearly identified 
three distinct phases in the (p, T)-diagram. In addition 
to the two known ones there is a crystalline phase with 
broken chiral and translational symmetry. These three 
phases are all separated by 2nd order lines which meet 
in one point. 

We first explain how we have obtained the boundary 
of the chirally restored phase. A 2nd order line separates 
the chirally restored phase from either the massive or the 
crystal phase, depending on p. Since the scalar potential 



vanishes continuously across this line, the line itself can 
be determined rigorously by using perturbation theory. 
We start from Eq. (4.1), add counter terms to the lin- 
early and quadratically divergent pieces of Eq. (4.4) and 
take the limit L — > oo, replacing the sum by an integral. 
Setting E(k) = k + A(fc) and linearizing in A(k) yields 



N 



2 

i2 



dfc 

2^ 



In 



(1 



-/3(fc-/x) 



)(1 



3 

Jo 



dk , , , . sinh 8k 
A( fc ) ■ 

27r cosh j3k + cosh (3p 



(5.1) 



Let us now assume that translational symmetry is also 
broken spontaneously and that — just like in Sect. Ill 
at zero temperature — only 5±i ^ 0. Using the results 
for ADPT of Sect. II, we find (with q = qi = n/a) 

\Si\ 2 



A(k) = q-k + sgn(fc - q)y/(k - q) 2 + |Si| 2 + 



2{k + q) 



A(k) = 



\Sl\ 



for k G 



1 1 

+ 



q 3q 
2' 2 



k — q k + q / 

for < k < q/2 ,k > 3q/2 . (5.2) 



We can actually expand the square root in \S\ \ 2 provided 
we are careful about the treatment of the resulting pole in 
the integrand. One can easily convince oneself that one 
has to take the Cauchy principal value part (this being 
the only remnant of ADPT). Dropping those terms which 
do not depend on Si then yields 



N 



1 



\Si\ 2 



£ dfc ( 



ln(2A) 
k 



(5.3) 



sinh/3fc 



k 2 — q 2 J cosh 3k + cosh Bp 



The phase boundary for the 2nd order transition under 
the constraint that q (or a) has a given value is obtained 
from the condition that the r.h.s. of Eq. (5.3) vanishes. 
For each q this defines a curve in the (T, ^i)-plane. After 
a rescaling of variables, 

n = /3k , pf = /3p , A' = 0A , q' = 0q , (5.4) 

we can solve explicitly the condition that the coefficient 
of |Si| 2 in Eq. (5.3) vanishes with the result 



clit (p',q')= lim 2A'exp $(//,<?', A') 



(5.5) 



where 

$(//,<z',A') = / 
Jo 



A' 



d/i 



sinliK 



(q 1 ) 2 — k 2 J cosh k + cosh p' ' 

(5.6) 
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FIG. 6: Perturbative determination of phase boundaries as 
envelope of a family of curves. Each individual curve corre- 
ponds to a certain wave number of the periodic potential. 




FIG. 7: Wavenumber characterizing the crystal period along 
the phase boundary, showing the presence of a "Lifschitz 
point" at [if 



If we plot T crit = l/(3 clit (n',q') against ^ crit = MV&rit 
for all possible values of q', the envelope of the resulting 
family of curves represents the sought for phase bound- 
ary. Incidentally, since we allow q to vary freely, we can 
also consider the limit q — > in which the period of 
the crystal becomes infinite. In this limit one recovers 
the results of the standard GN model solution which are 
usually derived by assuming So =/= 0. Thus our method 
enables us to deal with the translationally broken and 
unbroken situation on the same footing. 

In Fig. 6, we display a few representative curves ob- 
tained from Eq. (5.6). In particular the curve ABC 
corresponds to q' = 0. The section AB coincides with 
the 2nd order line of the standard solution, B being the 
tricritical point [3] . This part of the phase boundary sur- 
vives if we relax the assumption of unbroken translational 
invariance. 

The true phase boundary beyond the tricritical point 
can be generated as envelope of the curves shown in Fig. 
(6). For each \i > \i t there is one particular curve (la- 
belled by q) which touches the envelope at this ^-value. 
This allows us to define a function q(n) which describes 
how the period of the crystal depends on the chemical po- 
tential along the phase boundary. This is of some interest 
since q = ir/a is an order parameter for the breakdown of 
translational invariance. In Fig. (7) we show the depen- 
dence of the order parameter q on ji as one moves along 
the phase boundary. The solid line is the curve q = /i 
which is approached asymptotically by the full calcula- 
tion. At /i = ji t , the tricritical point of the old solution, 
we see a clear signal of a 2nd order phase transition with 
breakdown of translational invariance. The point B in 
Fig. (6) is therefore also a multicritical point in the re- 
vised phase diagram. It plays a somewhat different role 
though. In fact, it has all the characteristics of a "Lifs- 
chitz point" in condensed matter theory [19]. This type 
of multicritical point has been discussed in the context 
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FIG. 8: Revised phase diagram of the GN model in the (T, \x)- 
plane. The dashed line is the first order phase boundary from 
the previous, translationally invariant calculation. 

of magnetic or liquid crystals which exhibit both peri- 
odic and homogeneous ordered phases. At the Lifschitz 
point the wave vector of the periodic structure vanishes 
continuously, just like in Fig. (7). 

One can determine the envelope of the curves from Fig. 
(6) numerically and display it in the (T, /i) diagram. The 
result is drawn in Fig. (8) up to fx = 1.5. For larger values 
of the chemical potential, it is possible to derive the form 
of the phase boundary analytically (see Appendix) with 
the result 

T crit = ^ . (5.7) 

This curve does not intersect the T — axis, as expected 
on the basis of our T = results of Sect. III. 

From Fig. (5) at T = and Fig. (7) along the pertur- 
bative phase boundary one would expect another 2nd or- 
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FIG. 9: Revised phase diagram of the GN model in the (T, p)- 
plane. The dashed lines belong to the old phase diagram, 
where they enclose the mixed phase. This "droplet" region is 
superseded by the crystal phase featuring baryons. 
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FIG. 10: Entropy density versus chemical potential at T — 
0.1 showing the absence of latent heat, in contrast to the 
translationally invariant case. 



der line joining the Lifschitz point (T t , p t ) with the point 
(T = 0, /i = 2/tt). The scalar potential is constant to the 
left of this hypothetical line and periodic to the right, its 
period a diverging as one approaches the line from the 
right. Since the amplitude is non- vanishing on both sides, 
we can no longer invoke perturbation theory. We have 
determined this phase boundary numerically by comput- 
ing the value of the grand canonical potential under the 
constraint that there are 0,2,4 potential wells in the in- 
terval of length L. This is done along curves of constant 
T, varying fi in the relevant region. At the critical line, 
all three values of the potential should converge for suf- 
ficiently large L. Thereby, one can locate the value of /i 
where the instability with respect to breakdown of trans- 
lational invariance (kink-antikink formation) sets in. For 
technical reasons, this particular calculation was done in 
a much smaller interval (L~ 40 - 80) but the results are 
stable with respect to increasing L within the accuracy 
of our plot. The resulting critical line is also included 
in Fig. (8). The tricritical point of the old solution co- 
incides with the point where all three phase boundaries 
meet. The dashed curve is the first order line of the old 
solution [3] which has now become obsolete. It has only 
been included to highlight in which way the phase dia- 
gram changes. 

A perhaps physically more illuminating way of present- 
ing the phase diagram is given in Fig. (9) where we have 
transformed our phase boundaries into the (T, p)-plane. 
Here, the first order line of the old solution splits up into 
the two dashed lines which delimit the mixed phase re- 
gion (droplets of chirally restored matter in the chirally 
broken vacuum). This should be replaced now by the 
two solid lines going downward from the former tricrit- 
ical point and enclosing the crystal phase. In the new 
solution at T = the crystal phase is stable at all fi- 
nite densities. If one increases the temperature at fixed 



density, the outcome depends on the value of p. For 
p > p t — l^t/ft, when crossing the new phase boundary 
one goes directly into the chirally- and translationally 
restored phase in a second order transition. For p < p t , 
translational invariance is restored first. At some higher 
temperature the dynamical fermion mass vanishes and 
chiral symmetry gets restored as well. 

Our results suggest that the phase transition between 
crystal and massive Fermi gas is also a 2nd order tran- 
sition. Since the corresponding phase boundary has 
only been obtained numerically, we have computed var- 
ious thermodynamic observables along isotherms in the 
(T, ^i)-plane to check whether they are indeed continuous. 
In the old solution, quantities like p, s, u are discontinu- 
ous across the first order line. By contrast, we see no 
discontinuity in any of these quantities within our nu- 
merical accuracy. This is illustrated by way of example 
in Fig. (10) for the entropy density s as a function of p, 
(at T = 0.1). Whereas in the old solution the entropy 
exhibits a huge discontinuity and jumps to the chirally 
restored value of Eq. (4.6) at the first order phase transi- 
tion, the crystal solution interpolates smoothly between 
the massive and massless Fermi gas. It shows no sign 
whatsoever of a latent heat. The corresponding curves 
for p and u also favour a second order transition. Since 
here the variation near the phase transition is more rapid 
and it is difficult to do accurate computations very close 
to the phase boundary, these other observables are less 
well suited for judging the order of the phase transition 
and will not be shown here. 

Finally, we wish to point out that one can do more 
analytical work on the crystal phase in the asymptotic 
region of large p and large f3. As derived in the appendix, 
the relation between Ap,S\ and Ap,T in the crystal phase 
at large p, is the same as the one between m and T in 
the massive Fermi gas phase at p, = 0, see Fig. (11). We 
do not really understand the reason for this interesting 
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FIG. 11: Asymptotic behaviour of the amplitude of the pe- 
riodic scalar potential versus temperature, valid for large [i. 
The curve shown is identical to the dependence of the fermion 
mass m on T at /i = 0. 
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FIG. 12: Single particle energy spectrum for the NJL2 model 
with dynamical mass gap. 



kind of scaling behaviour. It may point to additional 
simplifying features which we have not yet fully exploited. 



VI. SUMMARY AND PERSPECTIVE 

Let us summarize our findings about the phase struc- 
ture of the GN model and try to understand the essential 
underlying physics. Together with our earlier results for 
the GN model with continuous chiral symmetry [7] (or 
equivalcntly the "two-dimensional Nambu-Jona-Lasinio 
(NJL2) model", we are now in a position to cover both 
variants of the model and compare their phase structure. 
So far, it has been assumed that the phase diagrams 
are the same for the discrete and continuous cases [18] 
and exhibit two phases: Massive and massless fermions. 
The only symmetry issue considered was the breakdown 
of chiral symmetry. The phase boundary in the {p, T)- 
diagram had lines of first and second order transitions as 
well as a tricritical point. Due to the first order transition 
a mixed phase could also appear. 

The central result of our investigation is the fact that 
both of these model theories can break translational in- 
variance as well. Since this is a continuous symmetry, 
the large N limit is of course instrumental for the phase 
structure [2, 6]. We will now argue that the dynamical 
mechanism behind the breakdown of chiral symmetry in 
the vacuum and translational symmetry in baryonic mat- 
ter is actually the same. 

Our point can most clearly be made at T = 0. Let us 
first recall the structure of baryonic matter in the NJL2 
model. The condensate is best pictured as a helix in 
the space spanned by (ifitp), (ipij 5 ^), and x (the "chiral 
spiral" [7]). Since the chiral winding number is equal to 
baryon number, each full turn of the helix correponds to 
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FIG. 13: Schematic plot illustrating the filling of single par- 
ticle states in the NJL2 model for vacuum (/j, = 0), matter 
(fi > 0) and antimatter (pi < 0), respectively. Dark dashed 
regions: filled orbits, light dashed: empty orbits, white: en- 
ergy gaps. 



one baryon. 

From the point of view of Hartree-Fock theory, it is 
instructive to take a look at the single fermion spectrum. 
In the vacuum one has dynamical mass generation with 
the spectrum sketched in Fig. (12). Due to the peculiar 
properties of this model, the whole picture drawn in Fig. 
(12) just moves upward or downward by an amount /j, at 
finite chemical potential. In Fig. (13), we illustrate the 
filling of single particle orbits for vacuum, matter and 
antimatter (at T = 0), respectively. We see that there is 
always a gap of the same width "floating" at the Fermi 
surface. This enables the system to make optimal use of 
the level repulsion at the gap for lowering its energy. The 
important dynamics is gap formation at the Fermi sur- 
face, which reduces to the surface of the Dirac sea in the 
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FIG. 14: Same as Fig. 12, but for the discrete chiral GN 
model. The gaps are due to the periodic scalar potential. 

special case p = 0. The system breaks whatever symme- 
try it takes to generate the gap, either chiral symmetry in 
the vacuum or chiral and translational symmetry in mat- 
ter (with one unbroken combination of the two). This 
scenario has a lot in common with the Overhauser effect, 
originally discovered in non-relativistic Hartree-Fock sys- 
tems [12]. The nice feature about the present relativis- 
tic application is the intimate relationship between mass 
gap and band gap. Note that in the NJL2-model the sin- 
gle particle spectra for matter (p, > 0) and antimatter 
{ji < 0) are different, although the energy densities agree 
(for the same This follows from the fact that under 
charge conjugation 

iW = 7VW, (6-1) 

the scalar potential is invariant but the pseudoscalar one 
changes sign. 

Now consider the discrete chiral GN model at T = 
0. Here the potential is purely scalar and identical for 
matter and antimatter. The single particle spectrum is 
sketched in Fig. (14) for the high density case (S±i ^ 
only), showing the appearance of a symmetric pair of 
gaps. (For lower densities, more gaps would appear, but 
they are not relevant for the argument.) The filling of 
the single particle orbits is indicated in Fig. (15) which 
now replaces Fig. (13). The common theme is evidently 
gap formation right at the Fermi surface. 

Turning to non-zero temperature the phase diagrams 
of the two GN models are very different, both from each 
other and from what had previously been thought. The 
continuous chiral model possesses two phases, the "chi- 
ral spiral" and the chirally symmetric massless Fermi gas, 
separated by a horizontal 2nd order line T = T C [2] . Tem- 
perature affects only the radius of the helix, chemical 
potential affects only the helix angle, and the two act 
independently. By contrast the discrete chiral model ex- 
hibits three distinct phases: Massless and massive Fermi 
gas and kink-antikink crystal. We have not been able to 
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FIG. 15: Same as Fig. 13, but for the discrete chiral GN 
model. As in Fig. 13, there is always a gap at the Fermi 
surface. 



compute the full phase diagram analytically in this case 
but had to resort to a numerical solution of the Dirac- 
Hartree-Fock equations. The phase boundary of the chi- 
rally restored phase and the grand canonical potential 
at large p, (3 can be understood perturbatively. We find 
that the tricritical point of the old solution is replaced 
by a Lifschitz point separating the homogeneous from 
the periodic ordered phase. Due to thermal effects which 
smear out the Fermi surface, the massive Fermi gas phase 
survives at higher temperatures and small p. 

We have not found any first order transition in our 
calculations. As suggested by Fig. (9) the crystal phase 
takes over the role played by the mixed phase in the 
old solution, although the quantitative details are quite 
different. As explained in Ref. [7], droplets of chirally re- 
stored matter with extra fermions can be interpreted as 
"bag model" baryons with a mass of Mb = N/y/2. This 
"wrong" baryon mass leaves its traces in the old phase 
diagram of the GN model, for instance in the value of 
p c where the first order phase transition takes place at 
T = (the droplets fill all space), or equivalently in the 
slope of the energy density £hf{p) at p = 0. In the re- 
vised phase diagram for the discrete chiral model, the cor- 
rect, lower baryon mass {Mb — 2N/-k) appears instead. 
The continuous chiral model has massless baryons and 
consequently the structure of the phase diagram changes 
qualitatively. Again, the behaviour of £hf{p) at p = 
(vanishing slope) reflects the presence of massless baryons 
and breakdown of translational invariancc is shifted to 
the point p = 0. This shows that we have achieved an 
overall picture of the phase structure of both GN mod- 
els which is now consistent with their different baryon 
spectra and therefore physically more reasonable. 
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APPENDIX: 

Asymptotics at large /x and (3 

At large [i and (3 one can use perturbation theory to 
determine S±\. Starting point is the grand canonical 
potential of Eq. (4.1). Single particle energies can be 
computed in ADPT, Eq. (2.13). Denoting Si = iSi by s 
and qi = ir/a by q to ease the notation, find 



N 
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2 r q ' 2 dk 
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2 f q dk, 
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q/2 2^ 
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(6.2) 



with 



s 2 

2 V A; + g ' fc — g 



1 1 

+ 



£ 2 = q-y/(k- q) 2 + s 2 + 



E% 



+ V(k-q) 2 + s 2 + 



2{k + q) 



2{k + q) 



(6.3) 



Asymptotically, we may use fi = —q, cf. Fig. (7). Then, 
integrals involving 



c -l3(E 2 -n) e -0(E 3 -n) 



(6.4) 



are exponentially suppressed and can be neglected. Inte- 
grals involving 



(6.5) 



are negligible in the first term of Eq. (6.2) (k — 0...q/2), 
while in the fourth term (k = 3q/2...A) one has to pull 
out the T — piece as follows, 

1 + e 0(E 1+ „) =e 0(E,+„) (l + c -P(E 1+ ^ ( g 6) 

Now, the T ^ part is negligible. Finally, in the 3rd 
term of Eq. (6.2), we again decompose 

! + e 0(E 3 +„) = e 0(E 3+ „) ^ + e -f3(E 3 +^ ; ( g ?) 



but have to keep both pieces. In this way we find 
9 ri 

N 



| T ^ln(l + c^)) 

P J q/2 2ir 

9 fill 2 Ah 



f 3c,/2 dk 
L ~" 

# + ^ + -ln(2A) S 2 



(6.8) 



Finally, the following approximations can be made: Eval- 
uate the integrals over E\ and E3 exactly and expand the 
result for small s up to second order, keeping the logarith- 
mic corrections. In the other integrals, drop the s 2 -term 
in E2,Es and extend the integration limits to (—00... q) 
and (q...oo), respectively. Keeping only the s-dependent 
terms in yields the simple final result 

W = ~ [l-21n(4s 9 )] 

_ JL f°dfcln(l + c-' 3 ^^ N ) . (6.9) 
n(3 J \ ) 

All of these approximations have been checked numeri- 
cally against the full result. 

Before proceeding, let us compare this result with the 
one which assumes unbroken translational invariance, 
but taken at T ^ 0, fi = 0, 



N 



m 

47T 



(1 -lnm 2 ) 



/ dfcln(l + c- /3Vm " +fc ") . (6.10) 

i"P Jo V > 

The structure is almost the same, since fx has been "eaten 
up" by q from the single particle energies in Eq. (6.9). 
Let us recall how to evaluate m(T) at fi — (here, the old 
solution was correct). Take the derivative of the grand 
canonical potential with respect to m and set it equal to 
zero, 



ln(m) 



dfc 







lo Vm 2 + k 2 \ e 0^+^ + l, 

(6.11) 

Rescaling the integration variable and setting rh = j3m 
then yields the solution 



(3 — m exp < 2 



dn 



1 



Vm 2 + k 2 e^™ 2+K2 + 1 



(6.12) 



If we evaluate the function on the r.h.s. and plot rh/(3 
against 1/(3, we recover the usual result for the tempera- 
ture dependent mass, vanishing at T c — e c /n. That the 
Euler constant appears in T c follows from the following 
useful relation [20], 



f 

Jo 



dq 



sja 2 + q 2 e V a2 +i 2 + 1 



= -W--\c + 0(a 2 ). 

(6.13) 
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If we now go back to Eq. (6.9) for the periodic case and 
minimize with respect to s, we find 



In(4sg) + 2 f 
Jo 



dk 



1 







(6.14) 

This can be solved explicitly for q as a function of s,(3. 
The result can of course only be trusted if q and j3 are 
large enough. We can reduce this problem to the preced- 
ing one (GN model at fi = 0) as follows: The relation 
between 4qs and AqT in the crystal model is given by the 
same "universal curve" as the one between m and T at 
\i = 0. 

We finally note the following consequences of our 



asymptotic analysis (/i > 0): The boundary line in the 
(s,fi) plane at large /i and T = is given by 



s = 



4fi 



(6.15) 



in agreement with Eq. (3.9). The boundary line in the 
(T, n) plane at s = has the asymptotic form 



(6.16) 
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